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NOTE ON STIRLING'S FORMULA 

By S. A. Cokey 

Stirling's formula,* written in the notation ordinarily employed for 
Taylor's formula, is 

/(« + x) =/(a) + *[f(a+x) +/(«)] -^? [f'(a + x) -/»] 
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Here B u B. x , • • • are Bernoulli's numbers, as follows : B l = 1/(5, B% =1/30, 
B 3 = 1/42, Bi = 1/30, i? 6 = 5/66, B e = 691/2730, B 1 = 7/6, etc. 
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be developed successively by Stirling's formula, simplifying each by the aid of 
the preceding, the following formula is readily obtained : 

/(« + *) =/(<■> + j^j {[/(« + «) +/<«>] 
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This formula (2) is admirably adapted for the numerical computation of 
functions for which the value of the first derivative for each of the successive 

* See Bertrand's Calcul differ entiel, p. 332. 
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and equidistant steps, a, a + x/m, a + 2x/m, a + Bx/m, • • ■, a + x, is 
known or may be readily found, or of functions whose successive derivatives 
increase rapidly in complexity with the order of the derivatives. The special 
value of the formula will be understood from the following considerations : — 

1. No odd-numbered derivatives, except the first derivative, enter into 
the formula and their computation is therefore not required, except in so far as 
may be necessary to determine the form of the even-numbered derivatives of 
higher orders. 

2. The values of the even-numbered derivatives need be found for only 
two values of the variable, viz., x = and x = x, the values of such deriva- 
tives for the intermediate steps having been entirely eliminated from the 
formula. 

3. The value of the general term may be made as small as desired by 
making m large enough, that is by finding the value of the first derivative for 
a sufficiently large number of intermediate points, corresponding to equidis- 
tant values of x. 

As (2 ) is but the sum of a finite number of series developed by Stirling's 
formula, it follows that (2) is convergent if each of the component series is 
convergent, and, consequently, that the same test which determines the con- 
vergence of each of the component series determines the convergence of (2). 

The following examples will serve to illustrate the usefulness of formula 
(2) for numerical computation. In each example the value of m and the 
highest value of n used are given, and the results obtained are correct to the 
last place of decimals. 

f 1 dx 

(1) ^ -=0.6931471805 (m = 10, n = 3), 

/ sin cc 

(2) / ~^-dx = 1.67476 (m = 10, n = 2), 
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(3) J tan x dx =0.346573 (m = 10, n=l), 
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(4) jf 1 (l^t^dx = 0.84829 (m = 10, n = 1). 
Hitbman, Iowa, Octobbu 1903. 



